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Abstract. We present an analysis of our extensive data-set on IRAS 08544-4431. It is the first object we discuss of 
our newly defined sample of stars, selected for their position in the 'RV Tauri' box in the IRAS [12] — [25] , [25] — [60] 
two-color diagram. Moreover, our selection criteria included an observed excess in the L-band, indicative of a dusty 
disc. The SED of IRAS 08544-4431 shows a broad IR excess starting already at H. Our optical photometric data 
reveal some evidence for deep and shallow minima in the light curve and a pulsation time-scale of around 100 days 
with a small amplitude (AV peak-to-peak = 0.17 mag). Our CORALIE radial velocity measurements show that 
IRAS 08544-4431 is a binary system with a period of 499 ± 3 days and a mass function of 0.02 Mg. Moreover, 
IRAS 08544-4431 is detected in both the CO (2-1) and (1-0) mm-wave emission lines. The triangular shape of the 
weak CO profile confirms that part of the circumstellar material is not freely expanding but resides probably in a 
dusty circumbinary disc. Our chemical abundance analysis of a high resolution spectrum of high S/N reveals that 
a depletion process has modified the photospheric abundances to a moderate extent ([Zn/Fe]=+0.4). All these 
findings confirm that the F-type IRAS 08544-4431 is another good example of a binary Post-AGB star surrounded 
by a dusty disc. The Ha P-Cygni profile shows ongoing mass-loss with a very high outflow velocity, the origin of 
which is not understood. The strength and velocity of the Ha-absorption are modulated with the orbital motion; 
the maxima of both quantities (~ 400 km s~\ 5 A respectively) occur at superior conjunction. 
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1. Introduction 

A new set of probable post-AGB stars was found in the 
IRAS Point Source Catalog by Lloyd Evans (1999 and 
in preparation). Their recognition follows the finding 
that the brightest of the RV Tauri stars contained in the 
General Catalogue of Variable Stars have IRAS photome- 
try which includes a reliable 60 /im measurement, so that 
they may be plotted in a [121 — [25] , [251 — [60] two-color 
diagram. iLlovd Evans! <)l985[) and iRaveendranl l)l989(l 
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remarked on the fact that the region containing these 
data lies in a well-defined and relatively thinly-populated 
part of the diagram, which corresponds to the greater 
relative brightness of these stars at 60 /im compared 
to the majority of the M-type variables which have 
extensive dust shells. This made feasible a search of the 
IRAS photometry for new examples of RV Tauri stars, 
selected by their mid-infrared colors rather than by their 
variability in blue light. A list of candidates was prepared 
by calculating the colors of all IRAS stars with reliable 
(quality 3) photometry at 60 /im and selecting those which 
were located in a rectangle in the IRAS two-color di- 
agram which contained most of the known RV Tauri stars. 

The defining rectangle is 

f [12] -[25] ee 1.56 + 2.5 log[F(25)/F(12)] = 1.0- 1.55 
\ [25] -[60] ee 1.88 + 2.5 log[F(60)/F(25)j =0.20- 1.0 
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Table 1. Basic parameters of IRAS 08544 (The coordi- 
nates are from the GSC). 



IRAS 08544 


Coordinates 


O2000 


08 56 14.2 




#2000 


-44 43 11.0 


Galactic 


1 


265.51 


coordinates 


b 


+0.39 



Over 1000 sources south of 20deg N declination 
were located on the Palomar, ESO and/or SRC Sky 
Survey prints. Those where a reasonably bright star was 
found were selected for visible-light spectroscopy and 
near infrared (JHKL) photometry. Initial spectroscopic 
observations made with the Unit Spectrograph and 
RPCS detector system at the Cassegrain focus of the 
1.9m Radcliffe reflector at the Sutherland observing site 
of the South African Astronomical Observatory covered 
the approximate range 3500 - 6700 A, with a resolution of 
7 A. This survey revealed stars of a wide range of spectral 
types with the excess at L i ndica tive of the presence of 
a dusty disk i Llovd Evansl Il997h . Those objects with 
spectral types, based primarily on the hydrogen lines in 
view of the many spectral peculiarities of RVTauri star 
spectra, of F, G or K were regarded as possible RVTauri 
stars and were selected for additional observations. These 
included UBVRI photometry, of those stars which were 
not too faint and did not suffer crowding problems, made 
at Sutherland as part of the SAAO service observing 
programme. 

Repeated near infrared and/or visible light photome- 
try showed that while some of the RV Tauri-like stars are 
variables of substantial am plitude, most are les s variable 
than typical RVTauri stars. ILTovd Evan! l)l999(l reported 
a correlation between amplitude and spectral type: the 
most variable stars have spectral types in the narrow range 
occupied by the bulk of the previously-known RVTauri 
stars, while those of lesser amplitude are mainly of ear- 
lier F subtype. This suggested that many IRAS detected 
RVTauri stars of the GCVS are simply those stars with 
dusty disks which happen to be passi ng through the inst a- 
bility strip of the Type II Cepheids l(Llovd Evanslll999h . 

We embarked on monitoring programs both in optical 
and near-IR photometry at SAAO and radial velocity at 
La Silla. These data were supplemented with single ob- 
servations to obtain high-resolution spectroscopy and CO 
mm-wave line emission. The aim is to characterize these 
evolved objects, compare them with GCVS RV Tauri stars 
and to discuss the relationship between the presence of 
such a disc and the eventual binarity of the central ob- 
ject. 

In this first paper devoted to this project, we report on 
the analysis of an extensive data-set obtained for one ob- 
ject: IRAS 08544-4431 (IRAS 08544 in what follows, see 
Table It is one of the stars in our sample with a 



smaller amplitude and is of spectral type F3, which was 
determined on the basis of a low resolution spectrum 
()Llovd Evanslll999lh The object was selected, despite the 
low quality flag (2) of the 60 /xm flux point, because of its 
LRS spectrum. In the literature this little studied object 
is referred to as a car bon star following the LRS IRAS 
spectral classification l|Kwok et al.lll997f) . The LRS spec- 
trum is, however, rather featureless and in this paper we 
argue that the spe ctrum is not that of a C-rich object 
l)01non et alJll98fl) and that the star is certainly not a 
carbon star. In following papers we will describe the over- 
all characteristics of the sample. In Sect. 2 we outline the 
observational dataset and the reduction procedures, next 
we discuss the Spectral Energy Distribution (SED). In 
Sect. 4 we show that IRAS 08544 is variable with a small 
amplitude and long timescale. There is some hint that 
the defining light curve of GCVS RVTauri stars is also 
present in IRAS 08544 , but a dedicated monitoring pro- 
gram is needed to study this in more detail. In Sect. 5 the 
radial velocity measurements are presented, proving that 
IRAS 08544 is a binary; we also used our CORALIE data 
to study the variation of the P Cygni profile of Hcv. In Sect. 
6 we present our CO line emission measurements and in 
Sect. 7 our chemical abundance analysis shows that deple- 
tion, which is a widespread phenomenon among RV Tauri 
stars, has slightly altered the photosphcric abundances of 
IRAS 08544 as well. In Sect. 8 we put our main conclusions 
together and end with a discussion. 



2. Observations 

2.1. Photometry 

Infrared observations on the JHK L system defined by 
the standard stars of ICarterl Jl990) were made with the 
Mark II Infrared Photometer on the 0.75 m reflector at 
Sutherland (SAAO Facilities Manual). L has the mean ef- 
fective wavelength of 3.5 /jm and not 3.8 /jm for which the 
same symbol is sometimes used. Standard errors for this 
bright source are 0.03 mag for J, H and K and 0.05 mag 
for L. Individual observations will be given along with 
those of other stars by Lloyd Evans (in preparation) . 

Photometric observations on the U BVRI system were 
made with the Modular Photometer on the 0.5 m telescope 
at Sutherland on 54 nights from JD 2450858 to 2451666, 
with repeated observations on several nights. Observations 
were continued w ith th e 0.75 m Automatic Photoelectric 
Telescope (APT) [|Mart 



tmcz e liD[2002) on 37 nights from 
JD 2451929 to 2452060. This concentrated run of obser- 
vations made with the APT proved invaluable in showing 
that the pulsation time-scale was of the order of 100 days. 
The observations were made with re ference to the Cousin s 
UBVRI standards in the E- regions l|Menzies et alJll989l . 
and reduced to the standard system using color terms de- 
rived from observations of standard stars of a wide range 
of color. There was no local standard star and hence the 
data include the increased random errors involved in a 
magnitude transfer through a large angle. The nominal 
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error of a single observation is typically 0.007 mag but 
may be substantially larger in U — B, especially on bright 
moonlit nights, to which has to be added the contribution, 
perhaps of order 0.01 mag, of the E-region transfer. 

2.2. Spectroscopy 

2.2.1. Low Resolution spectra. 

The first spectra were taken with the RPCS detector on 
the Unit Spectrograph at the Cassegrain focus of the 1.9 m 
Radcliffe reflector at the Sutherland observatory of the 
SAAO. Spectra taken over the range 3500-7000 A with a 
resolution of 7 A were classified against a set of MK stan- 
dard star spectra and a hydrogen-line type of F3 was ob- 
tained from the higher Balmer lines (since Ha at least is 
always in emission) in the same way as for a large sample 
of known RV Tau stars observed earlier with a photo- 
graphic detector on the same spectrograph. IRAS 08544 
is thus one of the stars, described as RV Tau-like on the 
basis of their spectral energy distribution and F-G spec- 
tral type, which lie outside and on the hot side of the 
instability strip (Lloyd Evans 1999). 




10 10 10 10 

Wavelength A (nm) 



Fig. 1. The SED of IRAS 08544 . For the optical measure- 
ments, the minimal and maximal fluxes are displayed. The 
photometric data are given in Table [21 The dereddened 
fluxes are shown using a total reddening of E(B — V)=1A5 
with the standard reddening law. Note the broad IR excess 
starting already at H. 



2.2.2. FEROS spectrum 

A high signal-to- noise and high resolution (A/<5A ~ 48000) 
spectrum was obtained with FEROS mounted on the 1.52 
m ESO telescope at La Silla (Chile). FEROS is a fiber 
fed echelle spectrograph which covers the complete op- 
tical region from 3700 to 8600 A. The spectrograph is 
equipped with a CCD camera incorporating a chip with 
2048 x 4096 pixels of 15 x 15 /im. The online reduction 
consists of background subtraction, order extraction, flat- 
fielding and wavelength calibration, using a Th-Ar spec- 
trum. We normalized the reduced spectra using smoothed 
splines on interactively determined continuum points. The 
FEROS spectrum was obtained on January 26th 1999. 
Sample spectra are shown in Figs. 1121 and 1131 

2.2.3. CORALIE data 

We have been monitoring IRAS 08544 with the spectro- 
graph CORALIE mounted o n the Swiss telesco pe at La 
Silla, ESO Chile. CORALIE l|Queloz et alJll999h is a high 
resolution (A/<5A = 50000) fiber fed echelle spectrograph 
which is constructed to perform high precision velocity 
measurements. Cross-correlation profiles are calculated 
on-line with a standard F mask. To improve the cross- 
correlation profiles we constructed a individual mask for 
IRAS 08544 tuned on the basis of the FEROS spectrum. 
The radial velocities are obtained by Gaussian fitting of 
these profiles. We started this monitoring in November 
1998 and have gathered 107 data points over a total time 
span of 1435 days. 



2.2.4. CO measurements 

CO measurements were performed with the 15 m Swedish- 
ESO submillimeter telescope (SEST) at La Silla. A dual 
channel SIS receiver was used to observe simultaneously 
at 115 GHz (the J = 1-0 line) and at 230 GHz (the J 
= 2-1 line). At 115 GHz the half power beam width and 
the main beam efficiency of the telescope are 45" and 
0.7 and at 230 GHz 23" and 0.5. Typical receiver tem- 
peratures are 110 K (115 GHz) and 150 K (230 GHz). 
Three Acousto Optical Spectrometers (AOS) were used 
as backends. The 230 GHz receiver was connected to a 
high resolution spectrometer (HRS) and a low resolution 
spectrometer (LR1), the 115 GHz receiver to a second low 
resolution spectrometer (LR2). Five exposures were taken 
in the period from October 1999 to August 2001. The 
system temperatures for the CO (2-1) observations varied 
between 250 K and 350 K. The intensity scales, reported 
in Fig. OH and Table 01 are given as main beam brightness 
temperature, Tmb = T^/tjmb, where T\ is the antenna 
temperature corrected for atmospheric attenuation using 
the chopper wheel method and tjmb is the main beam 
efficiency. 

Additionally we mapped the environment of 
IRAS 08544 with low resolution spectra of (J = 2-1) 
CO on July 28th 2000. A map spacing of 11" was used. 
Finally, also one low resolution (J = 1-0) HCN spectrum 
was taken. No line was detected to an r.m.s. level of 
0.003 in Tmb- For all observations we used dual beam 
switching, where the source is alternately placed in the 
signal and the reference beam, using a beam throw of 
about 12' . 
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Table 2. The photometric data used to construct the 
SED. For the optical and near-IR photometry, the mean 
magnitude, the standard variation (a) and the number of 
observations (N) is shown. The IRAS and MSX fluxes are 
given in Jansky. The IRAS quality label and for MSX the 
percent error is shown. 





<m> 


a 


N 


A 


flux 


error 


u 


11.85 


0.06 


103 


12 /im 


180.30 




B 


10.31 


0.05 


103 


25 /im 


158.80 




V 


9.13 


0.03 


103 


60 fim 


56.25 




Rc 


8.14 


0.03 


103 


100 


28.43 


L 


Io 


7.11 


0.03 


103 


Bl (4.29 fim) 


80.17 


9.7 


J 


5.65 


0.02 


2 


B2 (4.25 fim) 


91.45 


9.1 


H 


4.67 


0.02 


2 


A (8.28 Mm) 


163.09 


5 


K 


3.51 


0.02 


2 


C (12.13 /im) 


173.50 


3 


L 


1.59 


0.02 


2 


D (14.65 (im) 


156.70 


4 










E (21.34 Aim) 


171.57 
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This near-IR excess is generally not observed in single 
post-AGB stars with an expanding dust shell. It is, how- 
ever, not unique and several post-AGB stars also show 
a similar SED with a broad-IR excess. All these latter 
objects turned out to be binaries and for these objects 
the near-IR excess is interpreted as coming from hot dust 
which is stored in a dusty disc. The indire ct and direct 
eviden ce for this interpretation is given in IVan WinckeJ 
(2001) and references th erein. For several objects, the disc 
is resolved (e.g. AGHer. ljura, et al Jl200ft> . 

We conclude that the SED of IRAS 08544 shows a very 
broad IR excess around a weakly variable F-type central 
star, which indicates the presence of a circumstellar disc 
instead of a freely expanding outflow. Its IR colours are 
very similar to those of the bulk IRAS detected RV Tauri 
stars. 



3. Spectral Energy Distribution (SED) 

All photometric data used in constructing the SED are 
given in Table[3 The optical coordinates of the Guide Star 
Catalogue optical object is some 19" west, and 7" north 
of the IRAS position and on the edge of the 95% posi- 
tional uncertainty ellips as given in the IRAS point source 
catalog. The near-IR images obtained with the PICNIC 
camera llClas.jllflflflj) in J, H and K did not show other 
near-IR sources which might be connected with the IRAS 
source. Moreover, the MSX l|Egan et al.l fl9 99) position 
shows a correct association with the optical source so we 
can safely assume the IRAS fluxes are associated with the 
optical component as well. 

The spectral energy distribution is shown in Fig. ^ 
The full line represents the Kurucz model with the param- 
eters obtained from analysis of the FEROS spectrum (see 
Sect. 7). The total reddening otE{B-V)=lAh is found by 
minimizing the difference between the dereddened fluxes 
in the optical and in the J-band and the model atmo- 
sphere. A rough estimate of the interstellar component of 
the reddening can be obtained using the relation between 
the DIB strength at 5 780 and 5797 A and the reddening 
l|Krelowski et alj fl9991. There is a rather large scatter in 
this relation (correlation coefficient 0.8-0.89) but we ob- 
tain a E IS (B - V)=0.7±0.2. 

With [12]-[25] = 1.42 and [25]-[60]=0.75 IRAS 08544 
falls clearly into the 'RV Tauri' box of the colour-colour 
diagram which is shown in Fig. 3(c) of iLlovd Evans! 
ljl985h . Moreover, IRAS 08544 shows a very broad IR ex- 
cess, which starts already at H (1.65 jLim) and points 
to the presence of hot dust in the close environment 
of the star. Also this hot circumstellar dus t component 
is an observed feature for RV Tauri st ars l|Gehrzl Il972: 
iGehrz fc Nevl Il972t ILlovd Evansl 1198.1) . although com- 
pared to the RV Tauri stars studied in the last reference, 
the excess in H, K and L for IRAS 08544 is large. 



4. Pulsations 

IRAS 08544 is clearly variable, albeit with a small peak-to- 
peak amplitude of 0.17 mag in V and 0.23 mag in B. Since 
IRAS 08544 has the same IR properties as RV Tauri stars, 
we checked whether we for similar photospheric pulsation 
characteristics. The photometric data set did not enable 
us to recover a clear period. The long run of observations 
in a single season with the APT covers 131 days (see lower 
panel Fig. |2J); the shallow 'secondary minimum', which is 
a defining feature of the RV Tauri stars might be present 
(see lower panel Fig.EJ, but we lack a long enough contin- 
uous time series to label the light curve as due to a genuine 
RV Tauri pulsation. The pulsation time-scale seems to be 
typical for RV Tauri stars. The successive deep minima do 
not repeat well, however. The full data set has been an- 
alyzed using three period finding methods but no clear 
conclusions on pulsational timescales could be reached. 
We found two periods of nearly equal probability, 71.9 
and 89.9 days. Additionally, the observations were plot- 
ted against all periods at one day intervals from 67 to 
110 days as well as the double periods near 144 and 180 
days and periods of 50.4 and 216 days suggested weakly 
by the period-finding programs, and no other plausible pe- 
riods were found. The substantial scatter, with for exam- 
ple two exceptionally bright points (JD 2451608 and 9, see 
upper panel Fig. [5J somewhat remote in time from other 
observations, is observed also in genuine RV Tauri stars. 
The observational errors are not negligible compared to 
a V peak-to-peak amplitude of only 0.17 mag. The colors 
vary such that B — V and V — I are bluest when the star 
is bright; there is no significant phase difference and the 
amplitudes are only 0.02 to 0.03 mag. The U — B plot is 
dominated by intrinsic and observational scatter. 

IRAS 08544 is no genuine high amplitude RV Tauri 
star but is variable with a small amplitude on a similar 
timescale. 
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Fig. 2. The V-band photometric variations displayed by 
IRAS 08544 : upper panel : all photometry, lower panel 
: only the APT photometry. The amplitude is low (AV 
peak-to-peak = 0.17 mag) and the pulsation time-scale is 
near 100 days. 

5. Spectral monitoring 

The main purpose of our spectral monitoring is to ob- 
tain radial velocity measurements by cross-correlation of 
IRAS 08544 , hence the signal-to-noise of every spectrum 
is kept low. As an interesting byproduct, we can study 
the variability of large features, like the Ha-profile. In 
the spectra of RVTauri stars line deformation and even 
line splitting is observed. They are commonly ascribed 
to a passage of a shock moving through the atmosphere 
JOillet et ahlll 9891 iLebre Rr. Gilledll 99ltl7iiiTet et al.lll990l 
iFokinl l 1994). A cross-correlation profi le is a good trace r 
of this motion in the photosphere ijMaas et alJ 12002^1 . 
For IRAS 08544 deformations are sometimes present in 
the cross-correlation profiles (Fig. [SJ) , but they are small, 
indicating that the shock is weak. The Gaussian cross- 
correlation profiles obtained with a proper mask for 
IRAS 08544 result in precise radial velocity measurements 
(Fig. |2J) and symmetric profiles. 

5.1. Radial velocities 

The online reduction software for CORALIE provides 
cross-correlation profiles computed with a standard tem- 
plate of spectral type F. As we have a high signal-to-noise 
FEROS spectrum, we decided to use this spectrum to con- 
struct a proper template tuned for IRAS 08544 to increase 
the efficiency of the cross-correlation technique. The tem- 
plate spectrum consists of box-shaped emission lines with 
equal amplitudes. These lines were taken from the stan- 



-50 50 100 150 

v rad (km/s) 

Fig. 3. Comparison between the cross-correlation profiles 
obtained with a standard template (open circles) and a 
specific template for IRAS 08544 (full circles). 

dard template but only the lines, which are present in 
the FEROS spectrum and for which we could determine 
the FWHM, were selected for the IRAS 08544 template. 
We used this FWHM of the Gaussian fit as the width of 
the lines. In this way we reduced the original 3074 line 
mask to a mask of only 953 lines. In Fig. [3] we compare 
the cross-correlation profile, computed with the standard 
template, with the one computed with the IRAS 08544 
template. Qualitatively the IRAS 08544 profiles are more 
symmetric and the noise in the profiles is lower, resulting 
in nicer profiles. For both templates we fitted the cross- 
correlation profiles with Gaussians. The central veloci- 
ties of the IRAS 08544 template are shifted with -0.17 
kms _1 and a standard deviation of 0.24 km s with re- 
spect to the velocities of the standard template. For the 
IRAS 08544 profiles we obtained a mean FWHM of 9.31 
km s _1 and a mean depth of 0.11, for the standard profiles 
a mean FWHM of 7.01 km s _1 and a depth of 0.08. The 
IRAS 08544 template results, thus, in broader and deeper 
profiles than the standard template. 

We have accumulated 107 radial velocities, for which 
an orbital period of 499 days was found. Fig. 01 (upper 
panel) presents the radial velocities folded on the orbital 
period and the best least-square fit. The orbital elements 
are listed in Table 01 Our dataset covers 2.9 orbital cycles 
and the standard deviation of the O-C values for the whole 
data set is 1.06 km s _1 . For the eccentr icity we obtained 
a valu e of 0.14. Applying the criterion of lLucv fc Sweenevl 
ill97ll) . we find the eccentricity to be significant (p=6.9 
10 -3 %). Assuming an average inclination of 60 ° and a 
mass of 0.6 M for IRAS 08544 , we find for the unseen 
companion a mass of 0.3 M . The companion could be 
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Fig. 4. Upper panel : The heliocentric radial velocity data 
folded on the 499 day period. The full line is the radial 
velocity fit discussed in the text. Phase corresponds to 
periastron passage. Lower panel : the residuals for the ra- 
dial velocities for the 499 day period fit. 



Fig. 5. The ^-statistics of the PDM-analysis for the resid- 
uals of the radial velocities for a period of 499 days (see 
text). 5 bins and 4 covers were used in this analysis. 
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Table 3. The orbital elements of IRAS 08544. 





IRAS 08544 
a 


Period (days) 
a sin i ( AU) 
F(M) (Mg) 
K (km s" 1 ) 
e 

T (periastron) (JD245+) 

7 (km s _1 ) 

oo-c 

N 


499 3 

0.32 

0.02 

7.2 0.2 
0.14 0.02 
26 10 
1964 14 
62.5 0.1 
1.06 

107 



s 

I o 



1800 
JD (+2450000) 



a late-type main sequence star or a low-mass (helium) 
WD. A compact companion makes the evolutionary status 
of this object (with the non-zero eccentricity) even more 
puzzling. 

We have analyzed the residuals (Fig. 0] : lower panel 
and Fig. HO) f or periodicities . The ^-statistics of the PDM- 
method l|Stellingwerjll978j) are shown in Fig. [S] We sus- 
pect variability at a timescale of 100 and/or 30 days, but 
we cannot prove yet that this period is stable on a longer 
time scale. We interprete this variability as due to pulsa- 
tions. 



Fig. 6. The residuals from JD 2451605 to 2452000, show- 
ing a clear structure on a timescale shorter than the or- 
bital period. The time-scale of this low amplitude scatter 
is around 90 days which we interprete as due to the pul- 
sational amplitude in the radial velocity data. 

5.2. Ha-profile 

In all our CORALIE spectra Ha shows a strong P-Cygni 
profile of a fast stellar wind, which is highly unusual for 
an F-type object. The P-Cygni profile varies with time 
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(see Fig. (JJ and indicates that the mass outflow is vari- 
able or non-sphcrically symmetric. As a detailed model of 
the outflow is needed to derive accurately parameters such 
as the expansion velocity and the mass loss rate, we will 
only determine the order of the expansion velocity and de- 
scribe the observational characteristics of the profile. We 
took as the expansion velocity, the velocity of the intersec- 
tion of the continuum level and the tangent in the bending 
point in the absorption component of the P-Cygni profile, 
and found that it varies between 100 and 400 km s _1 ! The 
EW of the absorption component (integrated between the 
wavelength of the expansion velocity and that of the in- 
tersection of the Ha-line and the continuum level) varies 
between 1 and 5 A. Fig. [8] shows that both the expan- 
sion velocity and the EW of the absorption component 
of Ha are correlated with the orbital phase and reach a 
maximum between phase 0.13 and 0.26. This is around 
superior conjunction at phase 0.16 : i.e. when the primary 
is at its farthest position with respect to us and the sec- 
ondary at its nearest. Thus assuming the F-type star is 
the Ha emitter in the binary system , the mass outflow 
is not spherically symmetric but seems to reach a maxi- 
mum in the direction of the companion. The intensity of 
the emission peak varies from 1.3 to 1.8 times the con- 
tinuum level. As the Ha- emission intensity is very s ensi- 
tive to the mass loss rate IjKunasz \- Momsoij fliKv'S. the 
mass loss rate of IRAS 08544 seems to be variable during 
orbital motion. Finally, the variation of the EW of the 
emission component is smaller than that of the absorpion 
component and lies in the range of 0-2 A. In some spectra 
an extra small blue shifted emission or a slight recovery 
between two absorption components is observed (see up- 
per spectrum in Fig. |7|) . Just before publication two ad- 
ditional observations were inserted. They correspond to 
orbital phases 0.22 and 0.24 and do not show the high ex- 
pansion velocity and high absorption EW expected from 
the preceding cycles. However, their expansion velocity is 
still higher than those between phases 0.35 and 1.0. This 
could indicate that there is a orbital cycle to cycle varia- 
tion in the strength and velocity of the Ha-absorption. 

For the moment neither the origin nor the location in 
the system of the fast Ha-emitting gas is understood. 
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Fig. 7. The variable P-Cygni profile of Ha. The orbital 
phase is indicated for each spectrum. 
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Fig. 8. The expansion velocity (upper panel) and the EW 
of the absorption component of the P-Cygni profile of 
Ha (lower panel) versus the orbital phase. 



6. CO mm-wave line emission 

6.1. CO (J=l-0) and (J=2-l) emission lines 

The CO data for IRAS 08544 are shown in Fig. M For all 
spectra a baseline is subtracted. The random mean scatter 
in the spectra outside the line is of the order of 0.02 K. In 
Table^the central heliocentric velocities, the FWHM and 
the integrated intensities, obtained by Gaussian fitting, for 
the different lines are presented. 

The central velocity of the low resolution CO (2-1) line 
taken on 2/10/99 (Fig. [21 upper panel), which corresponds 
to an orbital phase of 0.94, is 61.3 km s -1 . However, our 
CORALIE data (see Fig. 0} point to a photospheric ve- 
locity of 70.6 kms _1 at this orbital phase. The variation 



Table 4. The central heliocentric velocity (kms^ 1 ) ,the 
FWHM (km s~ x ) and the integrated intensity (K km s _1 ), 
obtained by Gaussian fitting, for the different lines. The 
second quoted line is the central position of the map shown 
in Fig. [TUI 



line 


V c FWHM lco 


CO (2-1) LRS 2/10/99 
CO (2-1) LRS 28/07/00 
CO (2-1) HRS total 
CO (1-0) LRS total 


61.3 6.3 2.8 
61.3 6.4 3.1 
60.7 5.9 2.3 
61.7 5.4 0.9 



due to pulsation can not explain the difference between 
the orbital velocity and the velocity derived from the CO 
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line. This means that the CO emission is not following the 
orbital motion of IRAS 08544 . 

Comparing the two low resolution spectra of the CO 
(2-1) line, observed at 2/10/99 and 28/07/00 (Fig.® first 
and second panel respectively), the latter being the spec- 
trum of the central position of our map, we state that both 
have the same central velocities. This is valid for all our 
profiles. We can thus conclude that the central velocity of 
the CO lines is constant. It traces the motion of the system 
and not that of the primary star. However, a small differ- 
ence is found between the central velocity (60.7 km s _1 , 
taken from the averaged HRS CO (2-1) line) and the sys- 
tem velocity of 62.5 km s _1 derived from the CORALIE 
data. 

The CO emission of IRAS 08544 is weak compared 
to its 60 fim flux. For AGB-stars a correlation is ob- 
served between the integr ated intensity of th e CO (1-0) 
line and the 60 /xm flux l|Nvman et al.lfl992|) . Also sin- 
gle post-AG B stars like the 21 /im objects follow this 
correlation l|Van Winckel et aD Il999l) . Using this corre- 
lation for oxygen rich stars, the observed 60 /im flux of 
56 Jansky for IRAS 08544 corresponds to an integrated 
intensity of the CO (1-0) line of 5.5 K kms" 1 . This 
is a factor 6 higher than the observed CO (1-0) inte- 
grated intensity (lco =0.92 K kms -1 ). With respect to 
the 60 fim flux coming from the dust, the circumstellar 
environment of IRAS 08544 i s deficient in CO, which is 
common in RVTauri obj ects 1 Alcolea fc Buiarraballll99ll 



IVan Winckel et alJll99flli 

The triangular line shape of the CO lines is peculiar 
and differs from the CO line shape expected for a mass- 
losing (Post)-AGB star. For a spatially unresolved mass- 
losing (Post)-AGB star the CO line is rectangular for a 
low optical dep t h and parabolic for a high optical depth 
l|01ofsson et all Il993j) . For 89 Her, AC Her, HD 44179 
(the central star of the Red Rectangle), BM Gem and 



served dLikkel et al.l 


1 


991; Alcolea & Buiarrabal 1991; 


BuiarrabaletalJ 1988 




Jura et al. 


19951 I.Tura fc Kahanc 


1999|. All these objects show a weak and small CO emis- 



sion component with a FW HM in the range of .5-5 
kms _1 and for the latter four ljura fc Kahanel lfl999) ar- 
gue that these emission components are the signatures of 
long-lived reservoirs of gravitationally bound gas. With a 
FWHM of 6 km s" 1 , IRAS 08544 falls slightly out of this 
range, but a higher inclination angle of the orbital plane 
can increase the FWHM of the CO emission line. 

We conclude that the CO line-profile of IRAS 08544 is 
also indicative for the presence of gravitationally bound 
gas in a circumbinary disc. Note that due to the misclas- 
sification of the LRS IRAS spectrum, this object was ob- 
served and d etected also in the CO surv ey of infrared car- 
bon stars bv iGroenew egen et al ] J2002|) .' We checked the 
profile of IGroenewegen et alJ l)2002|) . the velocity is the 
same but they used the half-width-zero-intensity to com- 
pute the outflow velocity for which they derive 8.2 km s . 
We checked their line-profile which was kindly given to us 



but since we obtained a high S/N high-resolution spec- 
trum, we used our detection in what follows. 

To obtain a rough lower bound for the mass con- 
tained in the CO envelope w e follow the trai n of 
thought of I.Tura fc Kahanel ljl999h . I.Tura et alJ l|l997h and 
iKahane et alJ l|l998h . As an approximation we assume 
that the gas is in circular orbits around the binary sys- 
tem for which we adopt a mass of 0.9 M Q (0.6 M© for 
the Post-AGB star + 0.3 M Q for the companion). We es- 
timate the inner radius, Ri n of the disc from the CO line 
FWHM. If Verb denotes the orbital velocity of the gas, 

R m = GM/V D 2 rb (1) 

With V or b = 3 km s _1 from our CO (2-1) high reso- 
lution line, Rin = . The luminosity of IRAS 08544 
lies in the range of 3000-7000 L Q . We adopt here a value 
of 5000 L Q . At a distance of 90 AU from the star with a 
luminosity of 5000 L Q a blackbody grain has a predicted 
temperature of 250 K. 

If the gas is optically thin we can crudely estimate the 
column density of CO molecules in the beam, N(CO) 

N{CO)= /f!f% f Tmbdv (2) 



47T 3 /X 2 /lI/ 2 j 

where u is the dipole moment of the CO molecule 
(0.1098 D; IChackerian fc Tippin"glll983|) and v is the fre- 
quency of the (2-1) transition (230.538 GHz). In this equa- 
tion several assumptions have been made : (1) the mi- 
crowave background is ignored, (2) the partition function 
scales as kT ex and (3) the Planck function and the cor- 
rection term for stimulated emission can be simplified by 
linear expansion of exp(±hv / kT ex ) . We find N(CO) ~ 
6 x 10 15 cm~ 2 . 

The total mass of CO molecules (Mco) m the beam is 
the product of N(CO), the mass of a CO molecule (mco) 
and the projected area of the beam on the sky Abeam- For 
a beam assumed to be Gaussian with half-power diameter 
Ohpbw-, 



Abe 



am — nu HPBW 



D 2 /(4/n2) 



(3) 



with the distance to the star. If we adopt a dis- 
tance of 1 kpc, we find M co = 3.6 x I0 25 kg. With the 
assumption of [CO/H] = 10 -4 , the total mass in the CO 
envelope is 0.2 M Q . If the gas is optically thick, this value 
is a lower bound. 



6.2. CO map 

We mapped the environment of IRAS 08544 in low resolu- 
tion CO(2-l). The CO map consists of 14 positions with a 
spacing of 11" and is presented in Fig. 1101 The clear maxi- 
mum observed in the central position of the map, which is 
also presented in the second panel of Fig. El indicates that 
the CO line is not interstellar. We find that an asymme- 
try is present in the map: while positions (-1,1) and (1,-1) 
show a clear CO emission, the intensity of the CO emis- 
sion is lower in position (1,1) and CO emission is absent 
in position (-1,-1). This could mean that the CO emission 
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is anisotropic and extended. However, this asymmetry can 
also be caused by pointing errors in the different positions 
of the map, which can also explain the signal present in 
the position (2,0). 

7. Chemical composition 

7.1. Atmospheric Parameters 

The abundances are calculated on t he basis of a LTE 
model atmospheres of iKuruczl l|l993() and the LTE pro- 
gram MOOG of Sneden. A model atmosphere is deter- 
mined by its effective temperature (T e ff), the gravity 
(log g) and the metalicity ([Fe/H]). We determined these 
parameters on the basis of spectroscopic criteria : the ef- 
fective temperature by forcing the abundances of Fe I- 
lines to be independent of the excitation potentials and 
the gravity by forcing ionization balance between the Fe 
I and Fe II lines. Moreover, the microturbulent velocity 
was estimated by forcing the abundance of the Fel-lines 
to be independent of the reduced equivalent width. We 
then used the model atmospheres with the correct overall 
metalicity as given by the Fe abundance. The oscillator 
strengths of the lines were taken from the critically com- 
piled database continuo usly updated in our institute. Fo r 
more details we refer to IVan Winckel fe Revnierd l|2000|) . 

We obtained T eff = 7250 K, logg = 1.5 and [Fe/H] 
= —0.5 for the model parameters. This temperature is in 
agreement with the spectral type F3 determined from low 
resolution spectra and confirms the position of IRAS 08544 
on the high temperature side of the instability strip. 

7.2. Results 

The abundance results obtained for IRAS 08544 are listed 
in Table [3] In Fig. ^2 the abundance of the elements are 
plotted versus their dust condensation temperature. The 
spectrum has a high S/N of 150 (measured around 5800 
A) and we used a total of 264 lines, restricting ourselves 
to small lines (EW < 150 mA). These lines have, like all 
the other lines in the spectrum, symmetric profiles (see 
Fig. 1121 and I13|l . The good consistency between the abun- 
dances derived from the different ionization levels for the 
elements Cr, Mn, Ni confirms the reliable model atmo- 
sphere parameters. 

With [Fc/H]=-0.3 and [Zn/Fc]=+0.4 IRAS 08544 dis- 
plays a weak signature for depletion of refractory ele- 
ments in its photosphere. In Fig. 1111 we can distinguish 
two groups of elements. For elements like C, N, S, Zn, 
Na, Ba, Cu there is no correlation between the abundance 
and the condensation temperature. For elements with a 
condensation temperature higher than about 1200 K the 
abundance goes down from —0.2 for [Cr,Mn,Ni/H] to —1.1 
for [Zr/H], as the condensation temperature increases. 
One element doesn't rea lly fit in this picture which is Si. 
lLodders fe Feglevl l|l988|) give a condensation temperature 
of 1529 K for Si. However in the notes they say that most 
Si condenses into MgSiOz and Mg2SiOi which has a con- 



Table 5. Chemical analysis of IRAS 08544 . For the solar 
iron abundance we used the iron meteoric abundance of 
7.51. For the solar C,N and O abundances we adop ted 
resp. 8.57, 7.99 an d 8.86 (C: iBiemont et all Il993l N: 
Iffibbert et allll99l[ Q: lBiemont et al|l99lh.~r the solar 
Mg and Si abundances we adopt 7.54 l|Holwegerl20ol: the 
other solar abundances are taken from lGrevesse fe SauvaJ 
lll998l). The dust cond ensation temperatures are from 
lLodders fc Feglevl l|l988l) and references therein. They are 
computed using a solar abundance mix at a pressure of 
10~ 4 atm. They do not ascribe a condensation tempera- 
ture to O as it is the most abundant element in rock and 
therefore a separate condensation temperature is meaning- 
less. For all ions, the number of lines, the mean equivalent 
width, the absolute abundance and the a of the line-to-line 
scatter are given. 



IRAS 08544 
T eff =7250 K 
log g=l. 5 
^ = 4.50 km s" 1 
[Fe/H] '=-0.5 
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densation temperature of 1340 K. This lower temperature 
would improve the position of Si in our interpretation of 
Fig. EH 

To verify the depletion pattern observed in Fig. ^] 
we also calculated the abundances on the basis of two 
other model atmospheres. The first model has the param- 
eters T eff = 7500 K, logg = 2.0 and [Fe/H] = -0.5 , 
the second the parameters T e ff = 7000 K, logg = 1.5 



10 



T. Maas et al.: IRAS 08544-4431 



and [Fe/H] = —0.5. For the first model there is a global 
upward shift with a mean of 0.16, for the second model 
a global downward shift with a mean of 0.13. However, 
the typical depletion pattern is conserved for both mod- 
els. Moreover, the abundance ratios [S/a] (a — Si, Ca 
and Ti) and [Zn/Fe] which trace the depletion of a pho- 
tosphere, are almost independent of the model used. For 
the first model we obtained [S/a]=0.45 and [Zn/Fe]=0.42, 
for the second [S/a]=0.54 and [Zn/Fe]=0.35 (compared to 
[S/a]=0.50 and [Zn/Fe]=0.43 for our preferred model). So 
we can conclude that the depletion is small but real. 

The small N-overabundance ([N/H]=+0.2) might be 
the consequence of a first dredge-up. As the C-abundance 
is solar and the s-process elements follow the depletion 
pattern, there is no evidence for dredge-up of material 
processed by nucleosynthesis on the AGB. 

For the s-process element Nd we could only derive 
an upper limit since even the strongest expected lines 
were not present in the spectra. The upper limit is there- 
fore derived assuming a fixed EW of 5 mA for the 
strongest line in a clear part of the spectrum. We used 
the line lists of the Vienna Atomic Li ne Database (Vald2 , 
|http: / /www.astro.univie.ac.at /vald7 , iKupka et alJ fl9991 
to look for these strongest optical lines. 

We also noted a P-Cygni profile for two Fel lines, at 
5576.089 A and at 6400.001 A. We also detected 9 emis- 
sion lines in the spectrum (7 of Fel, 1 of Cal and 1 of Til). 
All these lines have low excitation potentials. 

To summarize, IRAS 08544 is a star where the deple- 
tion process has slightly altered the photospheric abun- 
dances ([Zn/Fe] = +0.4) and for which there is no evidence 
for typical 3rd dredge- up chemical enrichement. 

8. Conclusions and discussion 

The main conclusions obtained from the analysis of our 
large data-set are : 

1) IRAS 08544 is not a carbon star as listed in the 
literature but an F-type post-AGB star with IR colours 
similar to the bulk of the IRAS detected RVTauri stars; 
2) The star shows variability with a small amplitude 
(AV=0.17) for which the time-scale (between 30-100 
days) is not repeating well; 

2) IRAS 08544 is a binary system with a period of 499 ± 
3 days and a significant non-zero eccentricity of 0.14 ± 
0.02; 

3) the SED shows a broad IR-excess, indicating beside the 
presence of cold circumstellar dust also the presence of a 
strong hot dust component indicative of a circumbinary 
disc; 

4) the weak CO emission line and its triangular shape also 
points to a circumbinary location with a low gas-to-dust 
ratio; 

5) The Ha-line shows a P-Cygni profile for which both 
the strength and outflow velocity are correlated with the 
orbital period with maxima at superior conjunction; 

6) the photosphere is moderately depleted in refractory 



elements ([Zn/Fe]=+0.4). 

IRAS 08544 was selected on the basis of its IRAS 
colours whic h fall in to the RV T a uri b ox as defined by 
ILlovd Evans! lll985h . iKwok et alJ lll997h classified IRAS 
sources on the basis of their IRAS low-resolution spec- 
tra. They report the detection of circumstellar 11 /mi SiC 
dust emission for IRAS 08544 and classify the star as a 
possible carbon star. Our opical data clearly identifies the 
star as an F-type object and also the oxygen rich nature of 
the photosphere derived in our chemical abundance anal- 
ysis does not agree with a post-car bon star spectrum . 
Examination of the LRS spectrum l)Qlnon et all Il986|) 
shows that the 18 /jm silicate band is present in the form 
of weak emission, while the structure near 10 /im may be 
understood as the 9.7 /J,m silicate band with central rever- 
sal and low contrast. The spectrum therefore shows self- 
absorptio n. An alternative c lassification on the original 
scheme of lOlnon et all 1 198fi^ would be 31, indicating sili- 
cate absorption. Silicate self-a bsorption is explicit ly taken 
into account in the system of ILlovd Evansl l(l990h . where 
the appropriate classification wou ld be 21BA. The equ iv- 
alent on the new KSPW system l|Kraemer et al.ll2002f) is 
3. S B or 4. SB. This rec lassification removes the anomaly. 

ILlovd Evansl l)l999l) determined an RVTauri instabil- 
ity strip based on the RVTauri stars listed in the GOVS. 
The amplitude of IRAS 08544 is small for an RV Tauri star 
and its spectral type places it outside and on the hot side 
of this instability strip. This suggest that the star has al- 
most left the classical instability strip towards the phase of 
the planetary nebulae. Despite its probably more evolved 
nature, the infrared colours of IRAS 08544 are very similar 
to the IRAS detected RVTauri stars. 

Moreover, IRAS 08544 is a binary with a period of 
499 days together with a significant eccentricity of 0.14 
which is difficult to explain in terms of a previous standard 
binary AGB-evolution. On the AGB, the orbit was too 
small to accommodate an M-star with the same luminos- 
ity and it is not clear how the star avoided spiral-in, lead- 
ing to a circular orbit of short period. Post-AGB binaries 
are not uncommon and other objects with similar orbits 
and large eccentricities were observed (e.g. IWaters et alJ 
ll997tlVan" Winckcl 2003) . Population synthesis codes have 
shown that systems of post-AGB stars with resulting peri- 
ods in the order of one hundred to a few hundred days can 
only be expected if an efficient tidally enhanced mass-loss 
is induced. This wind increases the mass-loss significantly 
prior to Roche-lobe overflow (RLOF) and decreases the 
envelope binding energy to avoid spiral-in to very short 
orbits. Following IPols et all (|2003) no eccentric systems 
with P < 3000 d ays are produced; this is at variance with 
earlier reports bv lKarakas et, alJ I 2000F) . IRAS 08544 is an- 
other binary with a rather short period and eccentric orbit 
in an apparent post-AGB evolutionary phase for which the 
mere existence is not very well understood in binary stellar 
evolution theory. 

It has been known for some time that a subset of 
those post-AGB binaries are extremely depleted (e.g. 
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IVan Winckel et al.lll995h . The depletion process implies 
that the star is coated with a gaslayer, from which the 
dustgrains were removed, probably by radiation pres- 
sure. Since dust formation is thought to occur in the 
circumstellar environment and not in the photosphere, 
the process involves re -accretion of circumstellar gas. 
iMathis k, Lamersl l)l992f) have discussed several processes 
in which separation of circum stellar dust and gas may op- 
erate and lWaters et all 1 1992h added that the most favor- 
able circumstances may occur if the circumstellar dust is 
trapped in a stable disc. For post-AGB stars the presence 
of a circum-system disc implies binarity and the observa- 
tional basis of this process was that all strongly depleted 
objects known at that time were found to be binaries with 
strong evidence for some of them that the dust wa s indeed 
in a flattened geometry ijVan Winckel et al"lll995[) . 

The binary prerequisite for the depletion process was 
disputed, however, after the r ecent finding that many 



| ||196,^ suh- 



RV Tauri stars are also depleted. |Pre 
divided the RV Tauri stars into three spectroscopic classes 
(RVA,RVB,RVC) : RVA have spectral type G and K and 
show strong metallic lines and have normal CN and CH 
bands. RVB have spectral type Fp(R), they are weak- 
lined but show enhanced CN and CH bands. RVC are also 
weak lined objects with spectral type Fp but with nor- 
mal CN and CH bands. There is now general agreement 
that the abundance anomalies are mainly due to a deple- 



cesses ( 


Giridhar et all Il994! 


|l99 


Si 120001 iGonzalez et all 


1997a b 


; 




1998 


lMaas_el^Lll2QQd). The 



strength of the abundance anomalies differs from star to 
star. However, all the stars of the RVB class are affected, 
while the RVA class contains stars which are less or not 
affected at all. The RVC stars, which have a low initial 
meta llicity [Fe/H1<-1.0, d o not show any sign of deple- 
tion l|Giridhar et all EoOO). In IRAS 08544, the depletion 



altered the photospheric chemical composition but only 
mildly so. This is well established by the significant ra- 
tios [Zn/Fe]=+0.4, [S/a]=+0.5. As we have evidence for 
a current outflow, no accretion of gas is going on. The frac- 
tionation process has thus occurred in the past. However, 
the atmosphere is only slightly changed in its composition 
and since IRAS 08544 shows a large and broad excess , 
this illustrates again that no correlation exists between the 
strength o f the current IR-exce ss and the severity of the 
depletion I Giridhar et alJ f2f)00). An explanation put for- 



ward for the moderate depletion of th e RVA RV Tauri stars 
is convection ijGonzalez et alll997a|) : a depletion process 
has occurred but mixing with deeper layers with normal 
abundances, has diminished the effects. Convection could 
as well have diluted the anomalies caused by the depletion 
process in IRAS 08544 , when IRAS 08544 was cooler and 
ha d a deeper convectiye enve lope. 

IVan Winckel et aD l)l999|) argued that many IRAS se- 
lected RV Tauri stars are suspected binary systems where 
the circumstellar material resides in a stable dusty disc. 
Such a disc was resolved in a prototypical RV Tauri 
star AC Her ij.Iura et alJl2000h which was found to be a 



binar y with a period of 1200 days ijVan Winckel et all 
Recently, radial velocity monitoring showed that 
two strongly depleted RV Tauri stars, R U Cen and SX Cen , 
are indeed members of binary systems l|Maas et alJl20 02). 
Moreover, the orbital period of SX Cen is around 600 days, 
which is similar to the period of the long term photometric 
variation in the light curve. This indicates that the long 
term photometric variation in SXCen is due to variable 
circumstellar reddening and that the photometric class of 
the RV Tauri stars has its origin in a difference in incli- 
nation angle : the RVb RV Tauri stars are seen near edge 
on and the variable circumstellar reddening is due to the 
circumbinary disc, which is in the line of sight. The RVa 
RV Tauri stars are seen fa ce on and no long term p hoto- 
metric variation is present l|Van Winckel et 2]ES99). The 
binary system IRAS 08544 clearly fits in this picture. The 
hot dust observed in IRAS 08544 was not ejected recently 
but in a previous phase in which it was stored close to the 
star. 

Another interesting analogue to IRAS 08544 is 89 Her. 
It is a semi regular variable (pulsation period = 62 days, 
lFernieil98ll) of spectral type F2. It is a binary system with 
an orbital period of 288 days and an eccentricity of 0.19 
l|Waters et al.lll993|) . The SED shows a broad IR-excess 
and the CO emission profile has a ve ry narrow central 
peak with a width of less th an 1 km s~ 1 l(Likkel et alll99lt 
lAlcolea fc Buiarraballll99lT) . Interferometric CO observa- 
tions of 89 Her reveal the p resence of two distinct shells 
l|Alcolea fc Buiarraball9 95). The outer shell lies at a dis- 
tance of about 5 x 10 16 cm from the star and corresponds 
to the mass ejection in the AGB-phase, which stopped 
more than 1000 years ago (using an expansion velocity 
of 7 km s _1 ). The inner shell has an outer radius of 10 16 
cm. The Ha-line shows also a P-Cygni profile with a n 
outflow velocity of 150 km s~ 1 l)Sargent fc Osmerl fl969) . 
The chemical an alysis of the photosphere performed by 
iLuck et alJ 0990) proves, howe ver, that no depletion pro- 
cess has taken place in 89 Her. IWaters et al ] (119931) pro- 
posed for 89 Her a model in which dust has been stored 
in a circumbinary disc while no dust is formed in the low 
density outflow. In the plane of the disc the ejected mate- 
rial is decelerated. In the polar regions there is a free out- 
flow. A similar model could be applicable to IRAS 08544 . 
However, our Ha monitoring indicates that the maximal 
outflow velocity is towards the companion, not the poles 
(in which case a constant outflow velocity would be ob- 
served). Emission lines are observed in the spectrum of 
89 Her, as they are in IRAS 08544 . These emission lines 
could be collisionally excited as the wind of the centra l 
star runs into the circumsystem disc l|Waters et aljfl9 93) . 

On the basis of the triangular shape and the weak- 
ness of the CO line IRAS 08544 is similar to other similar 
objects : AC Her, HD 44179 (the central star of the Red 
Rectangle) , BM Gem and EU And. In these objects (a part 
of) the CO emission is coming fro m an orbiting molecu- 
lar reservoir l|jura fc Kahaneil999]) . The objects BM Gem 
and HD 44179 show, beside a small emission component, 
also a broad weaker component, probably coming from a 
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bipolar outflow. As the CO disk component is extremely 
small, the two components are clearly distinguishable. It 
is tentatively suggested that the CO profile of IRAS 08544 
also consists of two components which are difficult to dis- 
tinguish. For IRAS 08544 the CO disc component is wider 
and is more difficult to separate from a possible wide com- 
ponent. 

Beside a weak and small CO emission component, 
there are also other observations for HD 44179 and the 
RV Tauri star AC Her which indicate that the physical and 
chemical conditions in a circumbinary disc are different 
from those observed in an outflow : ISO-SWS spectra re- 
veal a high abundance of oxygen-rich c rystalline silicates 
llWaters et al . 1998; iMolster et al . 1999) and a very strong 
millimeter continuum flux from large dust grain s is present 
l|van der Veen et al.lll99l IShenton et aljll995|) . IR spec- 
tra are needed to verify whether these features are also 
present in IRAS 08544 . The photospheres of these three 
objects are depleted, but with variable strength : [Fe/H] = 
-0.3 and [Zn/Fe] = +0.4 for IRAS 08544 , [Fe/Hl = 
-1.7 and [Zn/Fe] = +0.7 for AC Her l|Van Winckel et all 
Il998h and fFe/H1=-3 .3 and [Zn/Fe]=+2.7 for HD 44179 
(|Waelkens et al.l fl996h The conditions under which the 
depletion process (whether or not covered by mixing) is 
highly efficient, of average efficiency or not at all efficient, 
is still not clear and needs further investigation. 

We conclude that IRAS 08544 is an interesting post- 
AGB binary with many characteristics of RV Tauri stars 
and which evolved probably just outside the instability 
strip. More detailed study, especially at high angular res- 
olution of IRAS 08544 and other similar objects will help 
us to clarify the relation between the characteristics and 
evolution of the circumstellar environment, the binary na- 
ture of the object and the presence of the depletion pro- 
cess. Moreover, the evolutionary link between the different 
post-AGB binaries needs further investigation. 
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Fig. 9. CO data : upper panel : low resolution CO(2-l) 
spectrum taken on 2/10/99, second panel : low resolution 
CO(2-l) spectrum taken on 28/07/00, third panel : to- 
tal averaged of 5 spectra of the high resolution CO(2-l) 
line , bottom panel : total averaged spectrum of the low 
resolution CO(1-0) line. All velocities are heliocentric. 



Fig. 10. The map of the environment of IRAS 08544 in 
low resolution CO(2-l). 
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Fig. 11. The composition of IRAS 08544 against the dust 
condensation temperature of the chemical element. 
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Fig. 12. Sample spectra of IRAS 08544 and HR 1865. The 
spectra are velocity corrected. HR 1865 is a massive super- 
giant with the same atmospheric parameters. All lines in 
HR1865 are stronger than in IRAS 08544 except for the 
lines of S. 
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Fig. 13. Spectrum of IRAS 08544 around 5660 A. All lines 
in HR 1865 are stronger than in IRAS 08544 , but the dif- 
ference in Sc-lines is largest. 



